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ABSTRACT 



'"i^ ' Context. Keplerian accretion discs around massive black holes (MBHs) are gravitationally unstable beyond a few hundredths of 
I , a parsec, and they should collapse to form stars. It has indeed been shown recently that an accretion/star formation episode 
O ' took place a few million years ago in the Galactic center (GC). This raises the question of how the disc can survive in AGN 
| and quasars and continue to transport matter towards the black hole. 
Aims. We study the accretion/star formation process in quasars and AGN with one aim in mind: to show that a spectrum 
i ^ i ' similar to the observed one can be produced by the disc. 

Methods. We compute models of stationary accretion discs that are either continuous or clumpy. Continuous discs must be 
£t~) • maintained in a state of marginal stability so that the rate of star formation remains modest and the disc is not immediately 
^ ' destroyed. The disc then requires additional heating and additional transport of angular momentum. In clumpy discs, the 
momentum transport is provided by cloud interactions. 

Results. Non-viscous heating can be provided by stellar illumination, but in the case of continuous discs, even momentum 
f ' transport by supernovae is insufficient for sustaining a marginal state, except at the very periphery of the disc. In clumpy discs 

it is possible to account for the required accretion rate through interactions between clouds, but this model is unsatisfactory 
0^ ' because its parameters are tightly constrained without any physical justification. 

CD Conclusions. Finally one must appeal to non-stationary discs with intermittent accretion episodes like those that occurred in the 

o : GC ' but such a model is probably not applicable either to luminous hish redshift q ~ 01 to radio " loud quasars - 

^ , Key words, quasars - AGN - Galactic center - accretion disc - star formation 

X' 

& 1 1. "Historical" introduction This is why in a previous paper (Collin & Zahn 1999b, 

hereafter referred as CZ99), we developed an exploratory 

It is well-known that Keplerian accretion discs around model to show that in the pre sence of star formation, the 

massive black holes (MBHs) are self-gravitating at accretion disc would be heated by the stars and could thus 

distances larger than about 1000 Schwarzschild radii be maintained in a marginally stable state where the star 

(10 3 i?schw)- Consequently they become gravitationally f 0rma ti O n rate is strongly reduced and destruction of the 

unstable to fragmentation on a short timescale and can gaseous disc is avo ided. Besides, we predicted that the 

collapse to form stars (Paczynski 1978; Kolykhalov & stellar mass function wou i d be top-heavy due to continu- 

Sunyaev 1980; Lin & Pringle 1987; Shlosman, Frank & mg accret i on onto the forming stars, giving rise to massive 

Begelman 1989; Shlosman & Begelman 1989; Hure 1998; stars . We postulated that the momentum transport was 

Collin & Zahn 1999a, 1999b; Gammie 2001; Goodman caused by supemova explosions. We mentioned that this 

2003; Thompson et al. 2005, and several subsequent pa- model could be re i evant f or explaining the presence of hot 

pers discussing the Galactic center). Actually, Shlosman massive sta rs in the central parsec close to Sgr A*. Our 

& Begelman (1989) showed that unless they can be main- main concern in that paper was to show that a rate of star 

tained in a hot or highly turbulent state, the discs should f ormation consistent with the rate of supernovae explo- 

transform rapidly into stellar systems unable to fuel ac- sions could sustain the accretion rate. In this sense it was 

tive galactic nuclei (AGN). Since AGN discs live for long res trictive. We found a solution, but only in a relatively 

times (10 7 8 years), it was mandatory to understand how narrow range Q f radius (about a decade), and relatively 

they can avoid such collapse. far from the MBH ( at about lrj 6 R Schw )_ 

At that time there was no proof that stars, particu- 

Send offprint requests to: Suzy Collin (suzy.collin@obspm.fr) larly massive stars, exist very close to SMBs. Only recently 
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have starbursts been discovered in Seyfert nuclei on the 10 
parsec scale in two Seyfert nuclei (Davies et al. 2007). In 
1999, a few groups of young stars were observed in the 
nuclear regions of some nearby AGN, but not closer than 
a few tens of parsecs. At this distance, the disc is glob- 
ally unstable as its mass becomes comparable to that of 
the MBH, and gravitational torques produced by bars and 
spiral waves mix the gas copiously and can transport the 
mass. 

Concerning the Galactic center (GC), ideas have 
evolved considerably since the publication of CZ99, be- 
cause many young massive stars have been discovered in 
the central parsec (Eckart, Ott & Genzel 1999). At first, 
their presence was considered very puzzling, because it 
was thought that the strong tidal field of the MBH should 
inhibit star formation. Gerhard (2001) thus proposed that 
these stars were born farther out in a star cluster and later 
drifted close to the MBH by dynamical friction with back- 
ground stars. Several recent discoveries have made this 
scenario implausible. Nayakshin & Sunyaev (2005) showed 
that the maximum possible number of young low-mass 
stars (YSO), deduced from the comparison of Chandra 
observations of X-ray emission near Sgr A* with the Orion 
Nebula, is too small for the cluster inspiral model. In ad- 
dition, they showed that the initial mass function (IMF) 
of the stars near the GC is top-heavy. Therefore the 
"in situ accretion disc scenario" gained credit and gave 
rise to various discussions (Levin & Beloborodov 2003; 
Goodman 2003; Milosavljevic & Loeb 2004; Nayakshin 
2004; Nayakshin & Cuadra 2005; Nayakshin 2006 and 
2007; Levin 2007; Nayakshin et al. 2007). It seems clear 
now that star formation has taken place in the GC in an 
accretion disc and that the duration of this episode was 
rather short. 

However, one is still faced with a problem concerning 
AGN: no answer has been found so far on how star forma- 
tion can be limited to a small amount or even be avoided 
during long periods of time in the accretion discs of AGN. 
As mentioned above, one possibility is that a stationary 
disc is maintained in a marginally stable state for a long 
time. A second possibility is that the disc is not stationary, 
and like the GC, is prone to transient accretion episodes. 
We explore these two possibilities here. 

Actually we tackle the problem in a "reverse" way com- 
pared to what has been done so far. Because it was claimed 
(Sirko & Goodman 2003) that stellar heating is unable to 
maintain an accretion disc in a marginally stable state 
without producing a spectrum very different from that 
of a typical AGN, our basic assumption is that the disc 
emits a spectral energy distribution (SED) in agreement 
with the observations of AGN. We first assume that the 
disc is stationary and then check whether an acceptable 
amount of transport of angular momentum needs to be 
added to the standard viscous one to maintain the disc in 
a marginally stable state. We conclude that supernovae do 
not cause enough angular momentum transport in a large 
fraction of the disc. We also briefly examine whether the 
dumpiness that would develop naturally in a gravitation- 



ally unstable disc could provide the required "viscosity" 
through collisions of clumps, and we show that it could 
indeed be the case, but only if the size of the clumps is 
restricted to a narrow range of values. We conclude that 
one probably has to seek a solution in the framework of 
non stationary discs, and we look whether they can be 
accommodated with the observations of AGN. 

In the next section we discuss the lifetime of the 
episodes of star formation and accretion both in the GC 
and in luminous AGN. In Section 3 we recall the theo- 
retical basis of stationary accretion discs and of the star 
formation mechanism. Results concerning typical quasars 
and Seyfert nuclei are given in Section 4. In Section 5 
we develop a simple model for accretion in a stationary 
clumpy disc, and finally, the last section is devoted to some 
considerations about non stationary accretion discs. 

2. How long should the accretion / star formation 
episodes last? 

2.1. The case of the Galactic center and of low 
luminosity AGN 

For the GC, our present state of knowledge has been sum- 
marized in Morris & Nayakshin (2007). For the stars, it 
is now known that most of the young massive stars in 
the central parsec of the GC are located in a relatively 
thin disk (aspect ratio ^0.1) rotating clockwise, while 
a few others are found in a diffuse disk rotating counter 
clockwise and inclined with respect to the first one (Levin 
& Beloborodov 2003; Genzel et al. 2003; Paumard et al. 
2006). From the spectral study of these stars, Paumard et 
al. (2006) have deduced the composition of the clockwise 
and counter clockwise discs. Using the evolutionary tracks 
developed in Geneva with different IMF, they found that 
the IMF is flatter than the Salpeter mass function, con- 
firming a formation in a non-common medium and that 
the proportions of different stellar types in the two discs 
can be well fitted by an episode of star formation taking 
place about 5 million years ago. They excluded a continu- 
ous formation process lasting more than two million years 
on the basis of the large number of WR stars. This im- 
plies that the gaseous disc around Sgr A* in which star 
formation took place has disappeared or has been unable 
to form new stars for 5 million years. 

Concerning the gaseous component, the inner 3 par- 
secs around Sgr A* constitute a highly complex region. 
There are infalling streamers of atomic and ionized gas 
(the "minispiral" ) in the inner parsec, embedded into a 
ring of molecular gas (the "circum-nuclear disc" , CND) 
peaking at 1.5 pc, made up of about 30 dense clouds with 
masses on the order of 10 4 M Q each and core densities 
of 10 7 ~ 8 particulcs per cc (Christopher et al. 2005). The 
CND rotates at 110 km s _1 with some deviations, and its 
structure seems consistent with a warped disc in Keplerian 
rotation around the MBH. The minispiral seems to be the 
inner edge of the CND, and both structures are in strong 
interaction. 



S. Collin, J. -P. Zahn: Star formation in accretion discs 



3 



Thus the inner regions located at a distance smaller 
than one parsec must permanently receive gas from the 
molecular clouds. The inflow rate deduced from the atomic 
gas has been estimated as 3 10~ 2 M per year (Jackson et 
al. 1993). Such an enormous rate cannot reach the MBH, 
unless it is much more luminous than observed. On the 
other hand, it has been shown that the presence of an op- 
tically thick disc in the inner 0.03 parsec can be excluded 
by the absence of eclipses of the star S2 and by the ab- 
sence of IR radiation due to reprocessing of the stellar flux 
(Cuadra et al. 2003). So the major fraction of the inflow 
must be hiding somewhere, perhaps in very dense molec- 
ular clumps too small to be Jeans unstable, which circu- 
larize at a short distance from the MBH, waiting for the 
accumulation of enough material to trigger a new episode 
of accretion and star formation onto the MBH. Note that 
the reservoir constituted by the CND is emptied in at least 
3 10 6 years. 

It is difficult to deduce the present accretion rate of 
the MBH precisely (cf. Quataert et al. 1999; Quataert & 
Gruzinov 2000b; Bower et al. 2005), but it appears that 
this rate can be provided by the hot stars, with no need of 
an accretion disc at the parsec scale (Cuadra et al. 2005, 
2006). That these hot stars formed 5 million years ago 
and that the accretion disc at the parsec scale is presently 
absent prove that the disc has not lived continuously for 
more than one or two million years, and that no further 
star formation episode took place. But what was the min- 
imum duration of this episode? Had it exhausted all the 
gas of the disc immediately? Until now this question has 
not been resolved. On simple probability considerations, 
we can deduce that such star formation episodes should 
be relatively frequent in the life history of quiescent galax- 
ies. Indeed young massive stars whose kinematics is con- 
sistent with a circular stellar disc in Keplerian rotation 
around the supcrmassive black hole are now observed also 
in the inner parsec of the Andromeda nebula (Bender et 
al. 2005). This is problably the first of a series, and similar 
stellar rings will be discovered in the future in other quies- 
cent galaxies with the development of optical interferomet- 
ric techniques, showing that a large fraction of quiescent 
galaxies possesses such rings of massive stars, implying 
frequent starburst events very close to the MBH. 

On the other hand low luminosity AGN (LLAGN, 
10 41-42 erg s _1 ) are observed in a large proportion of 
galaxies (Ho, Filippenko & Sargent 1997), implying that 
accretion episodes occupy about 10% of the lifetime of qui- 
escent galaxies. Nayakshin (2004) suggests that LLAGN 
go through episodes of quiescence with outbursts of ac- 
cretion and star formation. It is thus tempting to identify 
these LLAGN accretion events with the star formation 
episodes that occurred in the GC and in M31. Since they 
took place a few million years ago, their total duration 
should be at least of one million years. And since many 
short episodes would contradict the well-defined age of the 
hot stars in the GC, we deduce that the episodes should 
have been continuous for one million years. It means that 
the discs should have survived during a much longer time 



than the formation time of the massive stars (about 10 4 
years, cf. Nayakshin 2006, and below), so the stars should 
have formed at a low rate, or at least they should not have 
prevented subsequent accretion. 

Two different kinds of models have been proposed to 
account for the accretion/star formation episode in Sgr 
A*. In the first one, a gravitationally unstable stationary 
disc is supplied continuously at a high and constant rate 
at least on a viscous timescale. The flow is self-regulated 
by star formation and leads to a starburst in the region lo- 
cated at 0.1 pc, while allowing a small proportion of gas to 
reach the inner regions and be accreted by the black hole 
(Thompson et al. 2005). This model is only marginally 
consistent with the maximum duration of the star forma- 
tion episode, which is shorter than the viscous time scale 
(> 3 10 6 years). 

In the second model, a rapid episode of accretion and 
star formation took place, leaving the disc exhausted of 
its gas, through star formation and/or through accretion 
by the MBH, after a time much shorter than the viscous 
time of the disc. Nayakshin (2006) proposes that the life- 
time of the gaseous disc in the GC could have been pro- 
longed by the fact that, once the first stars are formed, 
the heating they provide is enough to maintain the disc 
in a stable state during a shorter time than the viscous 
time. In this model, the accretion rate is not constant 
across the disc, and the disc is in dynamical and thermal, 
but not in viscous, equilibrium. If it would last longer, it 
would reach a viscous equilibrium where its radial struc- 
ture would be reorganized, and additional angular mo- 
mentum would be necessary in order to keep the disc in 
this stable, or marginally stable, state. This is compatible 
with the duration of the episode. However the initial con- 
ditions require a delicate fine tuning so all locations of the 
disc reach marginal stability at the same time. The stars 
must begin to form at a high rate immediately after the 
disc is settled. On the other hand, they must not form too 
rapidly during the formation of the disc itself, or the gas 
would be destroyed. The formation of a disc able to form 
stars is a relatively long process, as it should satisfy two 
requirements: 

- First, a geometrically thin disc. If the disc is for in- 
stance created by two molecular clouds colliding and los- 
ing their angular momentum, they should also lose their 
quasi-spherical structures, and the time it would take to 
become a thin disc in hydrostatic equilibrium with the 
MBH is on the order of their dimension divided by the 
sound velocity, i.e. about 10 6 years. 

- Second, a Keplerian rotation. This would also take about 
this much time. 

From this, the process is compatible with an episode last- 
ing about 10 6 years. 

It is interesting to note that the amount of stellar 
heating needed to quench subsequent star formation in 
Nayakshin's model is > 10 40 ergs s _1 , according to his 
computation where the stellar mass is assumed to be 10 -3 
the disc mass, and the stars are assumed to have a mass of 
0.1 M Q and to be accreting at their Eddington rate. The 
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stellar luminosity is about the same as the black hole ac- 
cretion luminosity, if at this epoch it was accreting at 0.01 
of its Eddington rate. Extending this result to LLAGN 
means that the accretion and the stellar luminosity should 
be similar, as it indeed seems to be observed. 

Note that during the present "quiescent" phase of Sgr 
A*, there could also be some small intermittent accre- 
tion episodes, as proved by the discovery by Revnivtsev 
et al. (2004) of a hard X-ray emission of Sgr B2. It was 
presumably emitted 300 to 400 years ago by Sgr A* and 
Compton-scattered on Sgr B2. This leads to concluding 
that Sgr A* was at this time a kind of "very low lumi- 
nosity AGN" (10 39 erg s _1 ). But in such episodes the gas 
should not transit through an accretion disc on the par- 
sec scale (or star formation would take place), and regions 
closer to the MBH should pick up the required material 
directly, most likely devoid of angular momentum. 

We conclude from this discussion that though there 
is presently no massive accretion disc around Sgr A*, 
there should have been one a few million years ago, and it 
should have lasted for about one million years if such star- 
burst episodes are identified with LLAGNs. Star formation 
could have taken place early and could have been inhib- 
ited later by stellar heating during the rest of the disc 
lifetime, where the luminosity was dominated by newly 
formed stars. These episodes are expected to occur many 
times during the lifetime of a galaxy. But the problem 
of the initial formation of the accretion disc is still open. 
Note that they should lead to a large population of stel- 
lar remnants (neutron stars and black holes) and of solar 
mass stars in the inner parsec around the MBH. If these 
episodes occur every 10 years, it would mean that there 
should be presently about 10 5 stellar remnants in the in- 
ner parsec of the GC, which is indeed compatible with the 
observational constraints (Deegan & Nayakshin 2007). 

2.2. The accretion / star formation scenario in 
luminous AGN 

The accretion process in quasars and in luminous Seyfert 
nuclei is probably quite different from that in the LLAGN: 
it is due to mergers and close interactions of the host 
galaxy with other galaxies. This raises two questions: is 
the accretion continuous or intermittent like in the GC? 
Does a disc exist still at the parsec scale? 

Concerning the first question, one can easily imagine 
that the accretion is continuous and even close to the 
Eddington rate during a time comparable to the interac- 
tion (10 8 years), which is also the growing time of MBHs. 
There are at least two arguments in favor of this possibil- 
ity: 

1. the presence of very massive BHs (M > 10 9 M Q ) in high 
redshift (z ~ 6) quasars, implying a growth time of a few 
hundred million years, which would not leave any room for 
intermittent accretion; but one could argue that in these 
objects the accretion process is different from AGN and 
less distant quasars; 



2. the simultaneous existence of extended radio lobes at 
the million parsec scale and of compact radio jets in the 
same direction at the parsec scale, in radio quasars and 
in FRI and FRII galaxiesQ. This requires a causal link 
between the two structures, implying a continuous activity 
during at least 10 7 years. 

Concerning the second question, we know that QSO 
and AGN accretion discs are gravitationally unstable be- 
yond 1000i? S chw (0.01 pc for a MBH of 10* M Q ). Since 
a mechanism is required to transfer mass from large dis- 
tances at a high rate, the accretion disc is certainly not 
limited to a dimension < 0.01 pc, unless low angular mo- 
mentum gas is provided in the form of a quasi spherical 
accretion flow. 

Several observations seem to prove the presence of 
dense gas in Keplerian motion up to much more than 
10 3 i? S chw Maser spots in NGC 4258 (Greenhill et al. 
1995) are observed at a few tenths of a parsec (corre- 
sponding to 10 5 i?schw) in perfect Keplerian rotation. The 
broad Balmer lines in Seyfert galaxies are emitted by a 
gas located at about 10 4 i?schw (the broad line region or 
BLR), and a part of the BLR is probably a flat medium 
illuminated by the central continuum and rotating with 
a nearly Keplerian velocity (in NGC 5548, Peterson & 
Wandel 2000; in Mrk 110, Kollatschny 2003), which could 
be made of clouds ejected from the disc (Murray & Chiang 
1997). If the near- infrared continuum in quasars and AGN 
is produced by the accretion disc, like the visible and UV 
continuum, as confirmed by the continuity in time lags 
and size of the emission region (in Fairall 9, for instance, 
see Clavel et al. 1989), its radius should be larger than 
10 3 i?schw (see below). All these results argue in favor of 
large accretion discs, though they do not imply that the 
discs are homogeneous or stationary. Note, however, that 
the rapid flux variations observed in quasars are not due 
to changes in the accretion rate in the outer disc, but to 
thermal or dynamical instabilities, or to sporadic changes 
in the accretion rate, taking place in the inner disc at a 
few tens of i?schw 

In the present paper, we thus assume that the gas ar- 
riving in the nucleus has a relatively large angular momen- 
tum and circularizes at a radius larger than 0.01 parsec. 
We thus first study a stationary disc, i.e. supplied perma- 
nently during at least 10 7 years at a rate corresponding 
to the luminosity of the MBH. In this case one should 
seek a mechanism able to avoid the destruction of the disc 
through star formation. It would occur if the disc is prone 
to only a modest star formation rate during this time, 
being for instance marginally stable or made of gravita- 
tionally stable clumps. Second, we examine the possibility 
that the disc, like the GC, is not fed permanently, and we 
discuss what consequences it would have on star formation 
and on the feeding of the black hole. 



The latter are assumed to be radio quasars seen edge on, 
according the the Unified Scheme, cf. Barthel 1989; this is con- 
firmed by the presence of several AGN signatures, like broad 
infrared lines. 
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3. Marginally stable stationary discs: theoretical 
considerations 

A geometrically thin accretion disc becomes gravitation- 
ally unstable when the Toomre parameter Q is smaller 
than unity (Toomre 1964; Goldreich & Lynden-Bell 1965). 
It is defined as 



VLc s 
vrGS 



2VT+C 

C ' 



(1) 



where is the angular velocity, £ the surface density 2pH, 
p being the density and H the scale height. In the second 
equation we have expressed Q in terms of £ = inGp/Q 2 , 
which measures the self-gravityQ. 

When Q is close to unity, the instability growth rate 
becomes very low, and, according for instance to Wang & 
Silk (1994) the time for star formation tf or m is 



^form — ^ 



(2) 



It is thus necessary for a stationary disc to be marginally 
stable (Q ~ 1) to avoid rapid transformation into stars or 
compact bodies. 

Goodman (2003) argues that, in order to feed AGN at 
rates close to the Eddington accretion rate, while main- 
taining the disc in a marginally stable state, the required 
level of energy feedback from star formation should be 
impossibly large. This argument was repeated in several 
papers and was used to dismiss a continuous episode of 
star formation in accretion discs. 

The essence of the problem lies in the amount of mo- 
mentum transport that the disc can sustain. The gas in a 
marginally stable disc must be transported at a relatively 
high velocity to maintain the required accretion rate, or 
it will be strongly heated. Goodman (2003) discusses sev- 
eral ways of providing the disc with an enhanced momen- 
tum transport, and concludes that they were all difficult 
to achieve, so the transport had to rely on standard tur- 
bulent viscosity. Consequently, Sirko & Goodman (2003) 
studied the case of discs with standard viscosity and show 
that any heating mechanism able to maintain the disc in 
a marginally stable state would lead to the a strong in- 
frared emission that contradicts the observed spectrum of 
a typical AGN. 

In CZ99, it was shown that supernovae can ensure an- 
gular momentum transport, while the newly formed stars 
provide modest complementary heating, both being suffi- 
cient for maintaining the disc in a marginally stable state. 
On the contrary, in the present paper we will try to ac- 
count for the observed spectrum, without imposing a rate 
of star formation inside the disc that is consistent with 
the rate of supernovae necessary to transport the angular 
momentum. 

Let us recall the problem here. The standard sta- 
tionary "a-disc" in its simplified vertically-averaged form 



(Shakura & Sunyaev 1973) is governed by a set of alge- 
braic equations allowing solution of the radial structure as 
a function of only 3 parameters: the mass of the MBH, M, 
the accretion rate, M, and a, the "viscosity parameter'!!. 

The equation accounting for the transport of angular 
momentum is given by 



Mf(R) = SirvpH = 6irapc s H 2 , 



(3) 



where H is the disc scale height, p the vertically aver- 
aged density, and c s the isothermal sound speed y/ P j p, 
P being the sum of the gas and radiation pressure; v is 
the kinematic viscosity given by the a-prescription, a be- 
ing a dimcnsionlcss number < unity; /(i?) accounts for 
the inner boundary conditions. It is equal to unity in the 
outer regions of interest for us, so we will ignore it in the 
following. 

The accretion rate also intervenes in the energy bal- 
ance between local viscous dissipation and radiative cool- 
ing: 



F rad - aT eS 



m 2 M 



(4) 



where fl is the Keplerian angular velocity, a the Stefan 
constant, and -F ra d the flux radiated by each face of the 
disc. 

The inner part of the accretion disc is assumed to be 
a standard a-disc, like in Sirko & Goodman. When Q 
reaches unity, it becomes marginally stable, and the pre- 
scription Q = 1 is imposed (hereafter we call it a "Q- 
disc"). Thus in the Q-disc there is the additional con- 
straint Q = 1 (corresponding to C = 2(1 + v / 2)=4.8) with 
respect to the a-disc. It leads to a given dependence of 
the density on the radius: p = \J\ + ( S7 2 /(27rG) oc R~ 3 . 
One must therefore suppress one of the equations in the 
set governing the a-disc. 

Here lies the fundamental difference between us and 
Sirko & Goodman. They chose to suppress the energy 
equation 31 assuming that "some feed-back mechanism 
supplies just enough additional energy to the disc to pre- 
vent Q falling below the minimum Q ~ 1." And they kept 
the momentum equation [3J In contrast, CZ99, following 
Collin & Hure (1999), kept the energy equation |H but 
suppressed Eq. [31 assuming that some additional momen- 
tum transport mechanism exists in the disc and maintains 
it in a marginal state. In CZ99, this momentum transport 
was identified with the effect of the expansion of super- 
novae remnants in the disc, and the contribution of the 
star heating to the energy equation was taken into ac- 
count. 

Why was the amount of extra heating so important 
in the Sirko & Goodman paper? Simply because, as they 
assume a turbulent viscosity in the Q-disc and moreover 



( is generally neglected with respect to unity, although it 
can be very large for Q smaller than unity. 



some authors, like Goodman, add another free parameter 
to distinguish between a viscosity based on gas pressure or on 
total pressure. We have adopted the total pressure prescription 
in our studies. 
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chose a low value for the a-parameter (a=0.01), the trans- 
port was very slow, and consequently the disc surface den- 
sity was high, requiring strong heating. In CZ99, the re- 
quired heating was much lower, as the transport mecha- 
nism was more efficient than a purely viscous one. Less 
gas - therefore less extra heating - was necessary. The 
stellar heating was locally stronger than the gravitational 
heating, but it was negligible compared to the whole disc 
emission. 

3.1. Mass transport 

Several causes can compete in transporting angular mo- 
mentum - supernovae, disc-driven winds or jets, stars cap- 
tured by the disc, cloud collisions, etc. We focus here on 
the transport of angular momentum related specifically 
to the star formation process, i.e. supernovae explosions. 
Since we assume that the disc is living for more than 
10 7 years, massive stars have indeed amply time to evolve. 
This mechanism was proposed first by Rozyczka et al. 
(1995). The idea is that, when a supernova explodes in 
the accretion disc and the shell expands, a fraction of the 
total momentum (P% t, about 10 43 g cm s _1 for a super- 
nova producing a kinetic energy of 10 51 ergs) delivered by 
the supernova is supplied to the disc. Their 2-dimensional 
simulations of a supernova exploding in a thin Keplerian 
accretion disc have shown that it leads to a redistribution 
of the angular momentum. This is because the leading 
hemisphere of the supernova shell has an excess of angular 
momentum compared to the disk, while the trailing hemi- 
sphere has a deficit of angular momentum. Consequently 
the leading hemisphere is driven towards the center and 
the trailing one towards the exterior. Rozyczka et al. made 
an analytical estimate of the angular momentum given 
to the disc, which agreed with their 2-D simulations. In 
CZ99 we questioned their estimation and proposed an- 
other one, which unfortunately lowers the efficiency of the 
process (see Appendix). Clearly, 3-dimension simulations 
are needed to settle the question; in the meantime we shall 
consider both prescriptions. 

In the following we divide the momentum transport 
into two parts, one provided by the standard turbulent 
viscosity (Eq.[3]), M v ; scous , and one given by the additional 
mechanism, M non viscous . 



3.2. Non viscous heating 

The shape of the SED in the near infrared band has always 
been a problem, as it is not well-fitted by the standard disc 
emission (see for instance Koratkar & Blaes 1999 for a re- 
view). It is clear that, besides viscous heating, another 
process contributes to the heating in the outer regions. 
This extra heating can have several causes: gravitational 
instability itself (Lodato & Bertin 2001), illumination by 
the central source (for instance if the disc is warped), to 
illumination by stars embedded in the disc or outside, etc. 
We do not distinguish between these different processes, 



because they can have the same parametrization in terms 
of an effective temperature as a function of the radius. If 
we hypothesize that it is due to stars, it will give us the 
maximum number of stars needed to achieve the heating. 
Note that the heating provided by stars outside the disc 
is a factor two smaller than if they were trapped within 
the disc. Even if the disc is not thick in the sense of the 
Rosseland mean, a surface density as low as 10 _3 g cm" 2 
(which is never reached in this problem) is indeed amply 
sufficient for absorbing the UV radiation coming from ex- 
ternal stars in a thin layer at the disc surface, and for 
reprocessing it as infrared radiation towards the interior 
of the disc (cf. Collin & Hure 1999 for a discussion of this 
effect). The only difference with stars embedded in the 
disc is that half of the reprocessed radiation is emitted 
towards the exterior. The number of stars is thus esti- 
mated simply by assuming that L non viscous (2L non v i scous 
if the stars are not embedded in the disc) is due to hot 
massive stars radiating 10 38 erg s _1 eactQ In Nayakshin 
(2006), the stellar heating was ascribed to low-mass stars 
(0.1 M Q ) in an Eddington- limited process of accretion; 
but this process lasts only 10 4 years, and thereafter, all 
stars are transformed into massive stars, so it is not valid 
for a stationary long-lasting disc. 

The non-viscous heating is taken simply as 



v non viscous 



non viscous 



ttR F, 



irui viscous: wh.6r6 -fnon viscous IS 



viscous 



8tt 



(5) 



3.3. Time scales 



Two time scales are of interest for this scenario. One is the 
viscous time, i v isc, equal to Mdi sc /M where M is given by 
Eq. [3l and Mdi SC is the disc mass, here computed simply 
as Mdi S c = 7ri? 2 S. The second is the time i accr it takes for 
the stars to acquire a high mass by accretion. 

Once a clump has collapsed (in a dynamical time of ~ 
l/fl) and a protostar is formed, it undergoes a mechanism 
of accretion and growth, such as proposed by Artymowicz, 
Lin & Wampler (1993) for stars trapped by the accretion 
disc. A star of mass m* at a radius R can accrete all the 
gas contained in its Roche lobe of radius 



V ~M ) 



1/3 



R. 



(6) 



The accretion rate is given by the shear velocity corre- 
sponding to the radius of the Roche lobe at each moment 
~ f2i?H and is thus 



in. 



ftR^pH. 



One finds the accretion rate, 
ni* 1 /m*\ 2 / 3 



M 



(7) 



(8) 



4 In CZ99, the question was treated in a slightly more de- 
tailed way as we considered the optically thin and the optically 
thick cases separately, and we assumed that the stars were em- 
bedded in the accretion disc. 
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where i accr = M/(pHR 2 Q), and finally at a time t 

3 



m* _ m init 
17 _ M 



t 



3f„ 



(9) 



where m^it is the initial mass of the star. 

However, we will not assume that the star can reach 
its Roche mass, which is very high (up to 1O 5 M0!), but 
that the accretion stops when it approaches the outflowing 
wind rate, which probably occurs when the mass is close 
to a few tens of M . The accretion rate is also limited 
by the Eddington rate. We call t accl the time it takes to 
acquire a mass of IOMq. 

3.4. Opacity 

Opacity plays a fundamental role in this problem, as it 
relates the mid-plane temperature T of the disc to the 
effective temperature T e g through the diffusion equation, 
which can be written as 



off 3t/8 + 1/2 + 1/4t' 
where r is the frequency averaged opacity 



r = pHn = 



(10) 



(11) 



and k the opacity coefficient in g _1 cm 2 . Since T deter- 
mines the radiation and gas pressures (cf. the appendix), 
and therefore the scale height of the disc, the whole disc 
structure depends on the opacity. In CZ99 we used a scal- 
ing law, assuming that K was a constant on the order of 
unity, which is only valid for T < 100K. In fact k varies 
considerably in the range of temperature of interest to us 
in this paper, i.e. around 10 3 K. 

Actually, one should make a distinction between the 
Rosseland and the Planck frequency-averaged opacities, 
the first being valid in the optically thick case, the second 
in the optically thin one. In Collin & Hure (1999), the dis- 
tinction was performed according to the computations of 
Hure et al. (1994). Here we decided not to use the Planck 
opacities, which would have complicated the problem un- 
necessarily, as there is almost always an optically thick 
solution in the regions of interest to us. We use interpola- 
tions of the mean opacities computed by Ferguson et al. 
(2005), which cover the whole range of temperature and 
density needed without any extrapolation. 

4. Marginally stable stationary discs: results 

We refer to CZ99 for discussions concerning the survival 
of the accretion disc to the perturbations caused by the 
stellar formation and the presence of hot massive stars in 
the disc or by the supernovae explosions. 

We computed 3 models. The first one corresponds to a 
luminous quasar (mass equal to 10 9 M Q , Eddington ratio 
equal to 0.5), the second to a luminous Seyfert nucleus 
(mass equal to 10 s M Q , Eddington ratio equal to 0.1), 



and the third to a weak Seyfert nucleus (mass equal to 
10 7 M Q , Eddington ratio equal to 0.01). The bolometric 
luminosities and the accretion rates are, respectively, 

— for the quasar: L— 0.75 10 47 erg s _1 , accretion rate = 
23 M Q yr-\ 

— for the luminous Seyfert: L= 1.5 10 45 erg s _1 , accretion 
rate = 0.46 M yr _1 , 

— for the weak Seyfert: L= 1.5 10 43 erg s _1 , accretion 
rate = 4.6 10~ 3 M Q yr" 1 . 

Since one aim of this paper is to show that the SED 
emitted by the disc heated by a non viscous process can 
agree with AGN observations, we first determine the rate 
of heating (i.e. aT 4 ff ) as a function of the radius, so as to 
obtain a spectrum in agreement with a typical AGN SED 
(the one given by Sirko & Goodman 2003). The theoretical 
SED is computed by assuming a local black body emission 
at the effective temperature T e g: 



vL v = 



167r 2 cos(i)hv 



RdR 



exp{hv/kT c ff) — 1 ' 



(12) 



where cos(i) is the inclination of the disc axis on the line of 
sight (in principle i < it/ A for a type 1 AGN according to 
the Unified Scheme) , and R ln and i? ou t are respectively the 
inner disc radius (here 3 i?schw) and the outer radius (here 
10 6 i?schw)- The factor 16-7T 2 accounts for the fact that the 
luminosity of an AGN is deduced from the flux observed at 
Earth, assuming that it is a point source radiating isotrop- 
ically. We take an inclination cos(i) = 0.75. We assume a 
Schwarzschild MBH with an efficiency i] = 0.057, and a 
viscosity coefficient a = 0.3 (this parameter is necessary 
not only for the a-disc, but also for the Q-disc to com- 
pute the viscous contribution to the heating and to the 
momentum transport). A good fit to the spectrum in the 
visible and infrared bands is obtained for T c s cx i?~ 7 with 
7=0.47, starting at a given radius, which we call the tran- 
sition radius, -Rtrans- In the quasar case, the radius where 
the disc becomes gravitationally unstable, -R ma rg, being 
small (300 Rsch), -Rtrans is chosen equal to i? ma rg- In the 
other cases, -Rtrans is smaller than i? ma rg- It is indeed log- 
ical that the heating due to stars begins to be smoothly 
efficient before i? m arg, since the stars can surround the 
disc or might have migrated towards the MBH (cf. CZ99). 
Note that the value of 7 is close to 0.5, which corresponds 
to a uniformly distributed source of heating when it dom- 
inate^. 

The method to compute the SED is oversimplified, but 
it has the advantage of being relatively correct for an op- 
tically thick medium in the visible and infrared ranges 
where the Q-disc is emitting. Indeed the critical wave- 
length A cr it corresponding to the beginning of the emission 
of the Q-disc is close to /ic/(4fcT e ff(i? marg )). T cr f(i? marg ) is 



For an illumination by the central source, 7= 0.5 can be due 
to a strong flaring or warping of the disc, or to a hot medium 
back-scattering the radiation of the central source (cf. Collin 
& Hure 1999). 
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Fig. 1. The computed SED for three typical models of a 
quasar and a Seyfert nucleus, corresponding to Figs. 2, 3, 
and 4, compared to a typical AGN SED. In all cases the 
internal and the external radii of the discs are equal to 3 
and 10 6 i?schw, respectively. The dotted line corresponds 
to the luminous Seyfert with a clumpy disc described in 
the second part of the paper. The grey zone shows the 
region where the SED is not computed correctly and is 
flatter than displayed here. 



~ 3000K. Thus A crit ~ 10000 Angstroms. This method is 
very bad for the UV emission produced by the a-disc, due 
to electron scattering opacity that distorts the spectrum 
towards high energy and makes it "flatter" . In the X-ray 
range, the emission does not stem from the accretion disc, 
but most probably by a hot medium lying above the disc, 
like a corona or distributed magnetic flares (Haardt & 
Maraschi 1991, 1993, and many subsequent papers). 

Figure Q] displays the three theoretical SED and the ob- 
served AGN spectrum. They agree reasonably well within 
the range of wavelengths allowed by the approximation of 
Eq. [T3 (i.e. above 3000 Angstroms), taking into account 
the differences of SED from one object to the other and 
the difficulty of defining a "typical" AGN spectrum. 

The results corresponding to the 3 typical cases con- 
sidered here are displayed in Figs. [H [31 and [3] and the 
model equations are detailed in the Appendix. These fig- 
ures show the dependence of several parameters on the 
radius: the optical thickness r, the effective temperature 
T e ff, the mid-plane temperature T, the opacity k, the den- 
sity p, the surface density £, the disc mass Mdi sc = 7r£i? 2 , 
the aspect ratio H/R, and the ratio (3 of the gas pressure 
P gas to the total pressure P gas + P r ad- The following figure 
(left of the fourth row) shows the value of Q in the a-disc, 
and the value of a g in the Q-disc, allowing a check that 
the junction between the two discs occurs when a c g (an 
"effective" viscosity parameter given by Eq. [3]) is equal 
to 0.3, the value of a in the a-disc, and when Q = 1 in 
the a-disc (except for the weak Seyfert case when we did 



not find any solution for the a-disc below Q — 3). Next 
are displayed the ratio of the radial to the sound velocity 
Vi-ad/cs, the ratio of the viscous to the total flux, the stellar 
luminosity L sta r (assuming that the additional heating is 
due to stars), the rate of supernovae explosions required to 
transport the angular momentum, and finally the viscous 
time t v i sc and the accretion time i a ccr in the Q-disc. 

Because of the additional heating, Q is increased in 
the a-disc that extends farther than in the absence of 
this heating. Of course this makes the implicit assumption 
that, when the disc is heated, it adjusts its structure to 
the new conditions, which means that its lifetime is longer 
than the viscous time (~ 10 6 years for the quasar and the 
luminous Seyfert, ~ 10 7 years for the weak Seyfert, see 
Figs. Id and©. 

An immediately striking result is that there are several 
solutions at the same radius in some regions, correspond- 
ing to an optically thick and an optically thin disc. This 
also occurs in the a-disc (cf. Fig. 3). Except in the first 
panel (left of the first row), we only show those solutions 
corresponding to an optically thick disc, since it is the only 
one to give the spectrum computed with Eq. 1121 Even the 
optically thick solution is not unique and there are some- 
times three optically thick solutions at the same radius, 
two with the temperature decreasing with the increasing 
radius being thermally stable, and the intermediate one 
thermally unstable. Note that the junction between the 
a-disc and the Q-disc corresponds to the stable upper so- 
lution. 

This phenomenon is well-known, for instance, in the 
case of dwarf novae. However, here it is not due to the 
"ionization instability" around 10,000K, but to dust sub- 
limation around 1,000K. The reason of such behavior is 
that the opacity k passes by a strong minimum (varying 
by many orders of magnitudes) around 1,000K with re- 
spect to its value at 100K, where it is dominated by dust, 
and 10,000K, where it is dominated by bound-bound and 
bound-free atomic processes. When the optical thickness 
r decreases, T (roughly equal to T 1 / 4 T c fr) also decreases. 
The radiation pressure (proportional to rT 4 ) and the gas 
pressure (proportional to T) also undergo a big decrease, 
leading to the decrease in the sound velocity, therefore in 
the scale height and in the surface density. And finally a e g 
(inversely proportional to c^E) can increase by 3 orders of 
magnitude! Thus k entirely governs the disc structure. 

Since one of the stable solutions disappears at a given 
radius, there should be a jump in the structure at this po- 
sition. However the jump cannot be abrupt, because there 
will be some radial heat transport through both radiation 
and turbulence, which are not taken into account in our 
computation. It is also highly probable that a supersonic 
turbulence develops and contributes to the transport of 
momentum in this region. Moreover, the Keplerian law is 
no longer valid close to the transition, as there is a strong 
radial variation in the pressure. It is clear that our solution 
in this highly perturbed region is non reliable. 

In the region where two stable solutions exist, the disc 
can also undergo changes in structure, passing through 
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Fig. 5. T c ff versus S for the model M = 10 8 M©, 
Eddington ratio 0.1, a — 0.3, and a dependency of T e g 
versus R: T e g oc R~~ 0A7 with i?t ra ns equal to 0.01 pc. The 
Eddington ratio varies from 10 -4 to 0.5 along the curves, 
which correspond from left to right to radii decreasing by 
a factor 2. To make the figure more legible, one out of two 
curves is in thick lines. The interruptions of the curves 
are due to the limitations imposed on the optical thick- 
ness and on the temperature. 



one solution to the other in a thermal time, which is much 
shorter than the viscous time. It would be interesting to 
see whether this thermal instability can lead to a "limit- 
cycle" like in dwarf novae. Figure \5\ shows the behavior of 
T e ff versus the surface density £ at different radii, for a 
disc with a BH mass of 10 8 M© and an Eddington ratio 
varying from 10~ 4 to 0.5. As expected from Figs. 3 and 
4, there are multiple solutions between 0.01 and 0.04 pc. 
But in contrast to the dwarf novae case, there are three 
values of T e g for one value of S, and not the opposite. 
The solutions with positive 5T e g/6H are stable and the 
one with negative <5T e ff/<5£ is unstable. The left stable 
solution corresponds to low optical thickness (but greater 
than unity), the right one to a large optical thickness. 
This situation cannot lead to a "limit-cycle" , as a stable 
solution exists always whatever the accretion rate (cf. Fig. 

In all cases, the study of this difficult problem is be- 
yond the scope of the present paper. It concerns the region 
located between 0.05 and 0.3 pc for the quasar, 0.01 and 
0.03 pc for the luminous Seyfert, and between 0.003 and 
0.01 pc for the weak Seyfert. Though we agree that it is 
unsatisfactory, we leave this part of the disc out of this 
discussion. 



Note first that the Q-disc is highly perturbed due to 
the opacity gap within about one decade of radius, as soon 
as the a-disc becomes gravitationally unstable. The main 
properties concerning only the non-perturbed regions are: 

- The mid-plane temperature, the density, and the surface 
density decrease with increasing radius. 

- In contrast, the mass of the Q-disc is about constant 
with radius and always negligible compared to that of the 
MBH; the disc is geometrically thin, with an aspect ratio 
on the order of 10 -3 to 10~ 2 . 

- cv c ff reaches high values, while the standard a cannot be 
greater than unity. The radial velocity varies accordingly, 
but it remains subsonic in the whole disc. 

- We can compare the viscous time with the time it takes 
for a star to acquire 10 M© (i acC r)- The viscous time is 
larger than the accretion time, owing to the large rate of 
mass transfer, meaning that as soon as stars are formed 
they can be considered as massive stars. 

Let us now discuss the most important properties re- 
lated to our problem, still only in the non-perturbed re- 
gions. 

4.1. Stellar heating 

The non viscous luminosity is always lower by at least two 
orders of magnitudes than the whole AGN luminosity. If 
we attribute this luminosity to a population of massive 
stars, it means that about 10 7 massive stars per decade 
of radius in the quasar case, 10 5 in the luminous Seyfert 
case, and 10 3 in the weak Seyfert case, are necessary at 
all radii. This is quite large, but not impossible if we con- 
sider that these stars can be located outside the disc. Such 
a population could correspond to a scaled version of the 
hundred massive stars in the GC, where the accretion rate 
during the accretion/star formation episode was probably 
respectively 5, 3, and 1 orders of magnitudes lower. Since 
these massive stars live for about 3 10 6 years, during the 
whole active lifetime of the object (say 10 8 years) there 
should have been 30 times more massive stars and there- 
fore stellar remnants. Thus there should be about 3 10 8 
M© in stellar remnants around the MBH in the quasar 
case, 3 10 6 M© in the luminous Seyfert, and 3 10 4 M© in 
the weak Seyfert case. Even accounting for twice as much 
stellar mass due to low-mass stars (recall that the IMF is 
top-heavy), it stays lower than the mass of the MBH itself 
in the quasar case, and it is completely negligible in the 
Seyfert case. Note, however, that these remnants should 
have migrated outside the accretion disc, or the gas would 
be dominated by the stellar mass. Finally we recall that 
the heating can be provided (and is certainly partly) by 
another source than newly- formed massive stars, such as 
irradiation by the central X-ray source. Finally, in the case 
of an irradiation by external stars, one should take the 
stellar contribution amounting to 2L star into account in 
the observed spectrum. This contribution should appear 
in the UV if it is emitted by massive stars. Though it is 
not negligible (it is equal to about 5% of the bolometric 
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luminosity), it would not be as detectable as the the disc 



emission. 



4.2. Rate of supernovae 

The rate of supernovae per decade in radius and per year 
corresponding to our prescription and to that of Rozyska 
et al. (1995) decreases with increasing radius, mainly be- 
cause H increases (more rapidly than the radius actu- 
ally), so the momentum given to the disc by one super- 
nova increases. The rates are generally about one order of 
magnitude higher with our prescription than with that of 
Rozyska et al.. 

The required rates of supernovae with the Rozyska et 
al. prescription at the beginning of the unperturbed region 
are for the quasar, the luminous Seyfert, and the weak 
Seyfert : 10, 0.1, and 10 -3 per year, respectively, (in the 
quasar case, it is the same with our prescription and that 
of Rozyska et al.). Since one supernova releases a mass on 
the order of 10 M Q in the disc as a stellar remnant or gas 
(of which a fraction stays inside the disc), it corresponds to 
an accretion rate higher than the assumed accretion 
rate! It decreases farther out, but only becomes accept- 
able at 10 pc from the MBH. Only in the weak Seyfert 
case is it comparable to the accretion rate, and is lower 
than the accretion rate (which is necessary for the gaseous 
disc to survive) at about 1 pc. Moreover, we recall that 
the Rozyska et al. prescription certainly underestimates 
the required rate, as explained in the previous section. 

The weak Seyfert case is interesting as a supply of mo- 
mentum transfer is necessary only for a radius larger than 
0.1 pc. At a smaller radius, the turbulent viscosity is high 
enough to provide the mass transfer in the Q-disc. This is 
due to the additional non viscous heating, which remains 
compatible with the observed spectrum and corresponds 
to a reasonable amount of massive stars. 

One can definitively conclude that supernovae are not 
able to provide the required momentum transport except 
in the outskirts of the Q-disc. In order for the inner accre- 
tion disc to survive star formation, it is necessary to seek 
another mechanism for momentum transport. Finally, the 
inner parts of the Q-disc are probably highly perturbed 
and out of equilibrium. 



5. Stationary clumpy discs 

5.1. Mass transport by cloud collisions: theoretical 
considerations 

Duschl & Britsch (2006) have shown that the gravitational 
instability itself may lead to turbulence and thus viscosity 
in the disc. The other way to get momentum transport, 
on which we concentrate here, is that a gravitationally un- 
stable disc should fragment into clouds whose typical size 
is equal to the disc scale height. Several people have stud- 
ied the "viscosity" resulting from their interactions, and 
different prescriptions have been proposed (Goldreich & 
Tremaine 1978; Stewart & Kaula 1980; Lin & Pringle 1987; 



Jog & Ostriker 1988 and 1989; Shlosman & Begelman 
1989; Gammie et al. 1991; Ozernoy et al. 1998; Kumar 
1999; Duschl et al. 2000; Vollmer & Duschl 2001; Vollmer 
& Beckert 2002, etc). 

We adopt a very simple approach here, assuming that 
the disc is made of identical clouds of dimensions Z c i that 
have a mean separation d c \. The basic idea is that, since 
the clouds have an orbital velocity depending on their dis- 
tance to the MBH, they are subject to a viscous torque 
able to transform orbital motions into turbulent motions. 
They lose a fraction of their kinetic energy through phys- 
ical collisions and gain kinetic energy by gravitational in- 
teractions with other clouds. Gravitational interactions 
occur when the clouds approach a few times their tidal 
radius c?h = i?(m c i/A/BH) 1 / 3 (see Eq. [B]where m* is re- 
placed by the mass of the cloud m c \). The change in veloc- 
ity during these encounters is equal to the shear velocity 
corresponding to this radius, fidn. 

When the dimension l c \ of the clouds is close to their 
mean distance d c \, physical encounters occur at the same 
frequency as gravitational encounters. For the steady 
state, this leads to a dispersion velocity: 



disp 



{Gm cl n) 2/3 . 



(13) 



When the size of the clouds is smaller than their mean 
distance, the mean time interval between two gravita- 
tional interactions i gra v is ~ fi (d c l/<fe) , while the 
mean time interval between two physical collisions i co ii is 
~ n^ 1 (d c \/l c \) 2 , so the balance between losses and gains 
writes as 



disp 



(Grocl Q)2/3^L „ (Gm cl fi) 2 / 3 (p) . (14) 

£grav \ ^cl / 

The existence of this dispersion velocity has two conse- 
quences. First, it corresponds to a horizontal random walk, 
associated with an effective kinematic viscosity (Goldreich 
& Tremaine 1978; Stewart & Kaula 1980) 



Vcft 



disp 



"grav 



1 



(15) 



where n grav = (^^j is the number of gravitational inter- 
actions during one rotation period. It leads to an accretion 
rate of 



M = 37rE av feff ~ 37rE E 



disp 



'grav 



1 



(16) 



where S av is the average surface density. 

Second, when it is an isotropic velocity dispersion 
(if the motions of the clouds are not restricted to a 2- 
dimensional structure, as is the case when the size of the 
clouds is smaller than the scale height of the disc), it in- 
tervenes in the hydrostatic equilibrium. One uses then the 
following prescription, which corresponds to replacing the 
thermal pressure in the hydrostatic equilibrium by the dy- 
namical pressure: 

Vdisp 1 



H~A- 



(17) 
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where A is a factor of order unity. The mass of the disc 
writes: 



disc 



ttR SL 



o 11 h i 

irR max 



H 



(18) 



We do not specify the heating mechanism, where we 
do not know what proportion is transformed into heat 
and radiated. As in the case of the continuous disc, we 
parametrize the heat flux by an effective temperature de- 
pending on the radius. Note that, if the coverage factor 
of the disc is smaller than unity, it should be taken into 
account as a factor multiplying the emitted flux in Eq. [T2] 

To these "macroscopic" equations, one should add 
the "microscopic" ones corresponding to the state of the 
clouds. EquationfTOlbetween the physical temperature and 
T e ff is still valid; but for the computation of the optical 
thickness, one must now take into account the number of 
clouds on a given line of sight, which depends on the cov- 
erage factor of the clouds (l c \/d c \) 2 (H/d c \). Equation ITTI 
thus becomes 



r = pH k max 



d c \ J dd ' 



(19) 



5.2. Mass transport by cloud collisions: applications 

It is beyond the scope of this paper to develop a fully con- 
sistent model, taking into account in particular the forma- 
tion and the destruction of the clouds. Our aim is simply 
to show that a typical accretion rate can be achieved by 
the viscosity due to cloud interactions, for a steady disc 
made of small spherical clouds of same masses and average 
distances. 

Case 1 

As for the continuous disc, we would like to fit the typ- 
ical spectrum of an AGN. In a first step we thus compute 
the coverage factor. It intervenes simply as a multiply- 
ing factor in Eq. [T2J so the value of 7 required to fit the 
observed spectrum stays the same as previously (0.47), 
and we also assume that the additional heating begins at 
-Rtrans- We adjust the viscosity coefficient a so that the a- 
disc becomes gravitationally unstable at a relatively small 
radius to better fit the observations (a=0.01). In the sec- 
ond step we impose the constant coverage factor found to 
fit the spectrum, as well as a constant ratio l c \/d c \. We 
solve the set of equations given in the appendix, and we 
compute the physical characteristics of the clouds (size, 
distance, temperature, density, etc) in order that the ac- 
cretion rate be equal to the imposed value. 

To lighten the paper we give only the results concern- 
ing the luminous Seyfert. The spectrum is shown in Fig. 
[TJ it corresponds to a coverage factor equal to 0.2. The 
structure of the disc is displayed in Fig.[5]for a ratio l c \/ d c \ 
equal to 0.2. 

Here again, as for the continuous disc, there are mul- 
tiple solutions between 0.1 and 0.4 pc, and consequently 
jumps in the radial structure due to the opacity gap. The 
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shown in Fig.Q] The parameters of the model are M = 10 s 
Mq, Eddington ratio 0.1, a = 0.01, and a dependency of 
T off versus R: T cS oc Rr QA1 , with i? tr ans equal to 0.5 10~ 3 
pc (Note that due a low value of a leads to a low value of 
R 



grav j 



The coverage factor is 0.2 and the ratio l c i/d c \ is 
equal to 0.2. The scarcity of points in some regions is due 
to the method used to sample the solutions. 



transition between the a-disc and the clumpy disc takes 
place at 0.03 pc, as shown by a cS = a = 0.01. The disc is 
quite massive, reaching a mass close to that of the MBH at 
10 pc, but it is very thin (H/R between 10 -4 and 10~ 3 ). 
As a consequence, Q is lower than unity, but the clumps 
are gravitationally stable, since their masses are relatively 
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Fig. 7. Same as the previous figure, but now not trying to 
fit the luminous Seyfert spectrum. The parameters of the 
model are M = 10 8 M Q , Eddington ratio 0.1, a = 0.01, 
and a dependency of T e g versus R: T c ff oc , with 

-Rtrans equal to 0.5 10~ 3 pc. The coverage factor is 0.2 and 
the ratio l c \/ d c \ is equal to 0.3. 



low (< a few Mq) and the Jeans temperature is smaller 
than the clump temperature. The viscous time is long, 
typically a few 10 7 years, but stays compatible with the 
disc being stationary during the lifetime of the AGN. 

Since the coverage factor is equal to 0.2, the disc ab- 
sorbs only 20% of the external or internal heating flux, 
the rest being radiated towards the observer. Owing to 
the additional heating required to fit the observed spec- 



trum, it means about 10% of the bolometric luminosity 
should reach the clumpy disc. This value is high, but it 
is not unrealistic if the disc is warped. This is indeed ob- 
served with the maser spots of NGC 4258, which actually 
constitutes a kind of "clumpy disc" . 

Case 2 

We can also consider another type of solution where 
the emission of the clumpy disc does not account for the in- 
frared band. The observed infrared spectrum should then 
be provided by another source than the accretion disc, 
such as hot dust in the torus. We assume a radius de- 
pendency T c g cx i?~ 3 / 4 , intermediate between a uniformly 
distributed source (7 = 1/2), and a central point source 
at a small height above the disc (7 = 1), and a=0.01 for 
comparison with the previous model. 

The results for the luminous Seyfert are displayed in 
Fig. [7] The covering factor is 0.2, and the ratio l c \/d c \ is 
equal to 0.3. The main difference from the previous model 
is the amount of heating, which now corresponds to only 
0.1% to 1% of the bolometric luminosity. The clumps are 
gravitationally unstable for a radius larger than 1 pc, so 
this clumpy disc cannot survive beyond this radius. 

Though they can give the required accretion rate, 
clumpy discs are not satisfactory because they are quite 
sensitive to the coverage factor and to the ratio of the size 
to the distance of the clumps, for which there is neither 
physical justification nor observational constraints. This 
is why we now discuss non-stationary discs. 

6. Non stationary discs in AGN 

A solution allowing accretion without too much star for- 
mation in AGN and quasars consists in an intermittent 
process comparable to what is operating in the GC, i.e. 
successive waves of accretion lasting less than the viscous 
time of the outer disc. There would thus be no need for 
another mechanism of transport than the viscous one. The 
accretion rate is given by Eq. [3] 

The problem with the initial state of the disc has al- 
ready been mentioned about the GC: Nayakshin (2006) 
assumes that that it settles as a marginally stable disc. 
We make the same assumption. Then we have to define 
the non viscous heating. There are two options: either it 
is absent and appears only when the stars begin to form 
(cf. Nayakshin 2006), or it is present at the onset of each 
episode, for instance if it is given by a previous generation 
of stars. Since massive stars last for about 10 7 years, this 
is quite plausible. We adopt here the second hypothesis, 
which allows us to impose the resulting spectrum like in 
the study of stationary disc. As a result, the radial struc- 
ture is the same as in the case of a stationary disc, with 
only the accretion rate and the viscous time different. 

Figure[8]compares i V isc and M in the stationary and in 
the non stationary cases for the 3 previous models. The ac- 
cretion rate in the non stationary case decreases outwards 
where it is much smaller than in the stationary case. But 
it is interesting to note that it is equal to the stationary 
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accretion rate at the transition radius between the a-disc 
and the Q-disc in the stationary case. At this position, 
the viscous time is about 10 6 years, while the non station- 
ary viscous time is longer in the outer regions of the disc. 
Thus, if the duration of the accretion episodes is near 10 6 
years, the inner a-disc will be in viscous equilibrium and 
stationary during the accretion episode, while the outer 
disc will not be out in equilibrium. This is an interest- 
ing configuration, since the emitted spectrum during this 
episode will be similar to the stationary one, i.e. close to 
the A GN spectrum. 

That we do not know the initial state of the disc and 
that there is presently no physical basis for the marginal 
state hypothesis are not too serious. Actually, the disc 
can be in any other initial state: it must only stay opti- 
cally thick to UV radiation, to be able to reprocess it into 
the infrared band and emit a spectrum in agreement with 
the observations. Even without any additional momentum 
transport, it could be easily gravitationally stable (Q > 1) 
owing to the stellar heating. In contrast, the flow can be 
gravitationally unstable when it starts up, and since stars 
are formed and begin to accrete after only 10 3-4 years, it 
becomes rapidly stabilized and able to reprocess the stel- 
lar radiation for the rest of the accretion episode. For each 
of these situations, we could give some numerical exam- 
ples, but we avoid this here to keep this paper from being 
too long. 

It is interesting to see the effect on non stationary 
discs of the thermal instability observed in a small re- 
gion of stationary discs. If the accretion rate increases in 
the outer annulus of this region, and if the initial repre- 
sentative point is on the left stable curve (cf. Fig. 5), the 
annulus will evolve through equilibrium states of increas- 
ing M into a viscous timescale. When M reaches the value 
corresponding to the maximum of this stable solution, the 
annulus will pass suddenly (into a thermal time) to the 
right stable equilibrium solution, triggering a strong per- 
turbation (non Keplerian motions, shock waves, etc) which 
will rapidly affect the whole unstable region. Inversely, if 
the accretion rate decreases, when it reaches the minimum 
values of the right stable curve of Fig. 5, the representa- 
tive point will jump suddenly to the left stable curves, so 
one would expect that the unstable region answers any 
perturbation of the accretion rate more rapidly than the 
rest of the disc. 

Intermittent accretion has implications for the growth 
of quasars: it means that the "duty cycle" represents only 
a small fraction of the cosmological time, as the ratio of 
the duration of the accretion episodes to the interval be- 
tween two episodes is less than, say 10%. For luminous 
(and therefore massive) high-redshift quasars, it is cer- 
tainly too small. We also mentioned that it would contra- 
dict the continuity observed in radio-loud quasars between 
the large radio structure and the jets on the parsec scale. 
For local Seyfert nuclei, it is compatible with the small 
proportion of Seyfert galaxies among all galaxies (1%). For 
a nucleus triggered by an interaction with another galaxy 
that lasts about 10 8 years, it means that it will only be ac- 



tive during 10' years. This implies that all galaxies should 
undergo at least 10 major interactions during their life to 
account for the observed number of Seyfert. 

In summary, it seems that intermittent accretion could 
easily account for both stellar formation and accretion 
onto the black holes in Seyfert galaxies, while at the same 
time preserving an AGN spectrum comparable to the ob- 
servations, but it does not solve the problem of high red- 
shift and radio-loud quasars. 

7. Conclusion 

There are similarities between the Galactic Center and 
luminous AGN: in both cases, accretion should take place 
and build a massive disc at about one parsec from the 
MBH. This disc is gravitationally unstable and prone to 
stellar formation. For the Galactic Center, the duration 
of the episode of accretion / star formation is too short 
for the disc to reach a stationary state. The suggestion of 
Nayakshin (2006) that the disc is maintained in a stable 
state by stellar heating can work if there is a mechanism 
forming a thin Keplerian disc in a time shorter than 10 6 
years. If the process is identical to that acting in LLAGN, 
this implies that it should be repeated every 10 7 years and 
that there should be a permanent a population of massive 
hot stars very near the MBH around quiescent MBHs and 
in LLAGN. This is indeed observed in M31 and in the GC, 
and it begins to be observed within ten parsecs in some 
Seyfert galaxies (Davies et al. 2007). 

Our purpose in this paper was to study the accretion 
rate / star formation process in quasars and AGN and to 
show that it can produce a spectrum similar to what is 
observed in AGN. We first considered the stationary case. 
To deal with the main problem for stationary discs, the 
transport of angular momentum, we focused on two mech- 
anisms related specifically to the gravitational instability 
leading to the formation of clumps and to star forma- 
tion, through supernovae explosions in a continuous disc 
or cloud collisions in a clumpy disc. 

We studied first the case of a continuous disc. Since 
the lifetime of the disc is longer than the formation and 
growth of the stars, the rate of star formation should be 
maintained within modest limits so that the gaseous disc 
is not immediately destroyed. It implies that during a time 
close to or longer than the viscous time, the disc should be 
maintained in a state of marginal stability (a "Q-disc"). 
However, it then requires additional heating and addi- 
tional transport of angular momentum. Starting from the 
observations, we studied the case of a disc accreting at a 
rate corresponding to the luminosity of the MBH. For the 
additional heating, what is provided by the newly formed 
massive stars can easily account for the spectral distribu- 
tion of AGN; and at the same time it contributes to main- 
taining the disc in a marginally stable state. However we 
have shown that the transport provided by the supernovae 
issued from the massive star formation is not sufficient for 
producing the required extra amount of momentum trans- 
port, except possibly in the weak Seyfert case. For quasars 
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and luminous Seyfert, it is sufficient only at distances ex- 
ceeding 10 pc. We have also shown that the Q-discs are 
perturbed by the existence of big jumps in their structure, 
due to the variations in the opacity with the temperature. 
It is probable that these jumps induce Shockwaves and 
supersonic turbulence. 

Then we examined a simple model of clumpy disc 
where the transport of angular momentum is provided by 
the interactions between the clouds. We showed that this 
transport is very efficient and can produce the required ac- 
cretion rate, provided that small clumps are formed and 
maintained inside the disc. These clumps are hot enough 
to be gravitationally stable. The observed spectrum can 
be accounted for by the model, provided that the disc re- 
ceives about 10% of the total luminosity of the object. The 
accretion rate can also be easily maintained if the disc is 
not irradiated as much, but in this case its infrared emis- 
sion is negligible. However, this model is unsatisfactory 
as the parameters (size and distance of the clumps) are 
tightly constrained without any physical justification. 

Finally we considered intermittent accretion similar to 
that in the GC. If the duration of the accretion episodes is 
of the order of 10 6 years, the outer disc is not stationary 
and not in viscous equilibrium, in contrast to the inner 
disc, whose accretion rate corresponds to the luminos- 
ity of the objects. As a result the spectrum emitted by 
the inner and the outer discs agree with the observations. 
This model could work for Seyfert galaxies, but neither 
for luminous high-redshift quasars (because the duty cy- 
cle would be too short), nor for radio- loud quasars (be- 
cause it does not ensure the required continuity between 
the large-scale and small-scale structures). 

A lot of work remains to be done. First, the structure 
of the "jumps" in radial structure due to the opacity gap 
should be examined in detail, since they are present in all 
kinds of discs, whether gravitationally stable or unstable, 
stationary or not, clumpy or not, and it is possible that 
the turbulence generated by these jumps has an impact 
on the angular momentum transport. Second, the clumpy 
disc deserves a thorough investigation, as it could be an 
interesting solution for the outer regions. And finally it 
would be necessary to study a genuine time-dependent 
model taking into account the evolution of star formation 
and its feed-back on the disc, as it could certainly repre- 
sent a solution for Seyfert nuclei. 
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Appendix A: Equations of the continuous disc 
structure 

We recall here the set of equations for solving the radial 
structure of the a-disc. Some of them have already been 
given in the text, but we prefer to rewrite them to get a 
fully consistent ensemble. We do not redefine the param- 
eters already mentioned in the text. 

The equation of hydrostatic equilibrium writes as 



c 2 = n 2 H(i + o 

with 

c s = — 

p 

pkT 



, and 



p — 



(A.l) 

(A.2) 
(A.3) 

(A.4) 



where p is the mean mass per particle, toh the proton 
mass. Then 



p 

1 rad — • 

2c 



where r is the Rosseland mean opacity: 
~2~' 



pHn 



(A.5) 



(A.6) 



where k is the opacity coefficient in g 1 cm 2 . 
The diffusion equation writes as 



aT 4 



3r/8 + 1/2 + l/4r' 



(A.7) 



The energy balance between the local viscous dissipa- 
tion and the radiative cooling requires 



3n 2 Mf(R) 
8^ ' 



(A. 



When there is an additional non viscous heating, this 
equation is simply replaced by T c g cx i?~ 7 , ensuring con- 
tinuity of T e ff at the transition radius. 

The cv-prescription defines the viscosity coefficient as 



(A.9) 



v being the kinematic viscosity. And finally the equation 
accounting for the transport of angular momentum is 




Mf(R) 
3tt 



= 2vpH. 



(A.10) 



When solved self-consistently with the appropriate 
functions for n(p,T) and p(p,T), these equations yield 
the a-disk structure. 

In the Q-disc, one has 



Q = 



n 2 



ttGS 2nGp 



Vi + C = i, 



(A.ll) 



which requires suppressing one of the equations in the 
previous set. As explained in the text, Sirko & Goodman 
(2003) suppress the energy-balance equation IA.81 while 
we suppress the equation for the moment transport IA.101 
However, as explained in Section 3.2, we take the viscous 
transport into account by allowing a fraction of M - which 
also contributes to the heating - be due to the viscous 
transport. The set of equations is solved by two different 
iterations for the a- and for the Q- discs. 

Since we are dealing with a gas between a few hun- 
dred and a few thousand degrees in the Q-disc of interest 
for us, we take fj, = 1, a value intermediate between 2.4, 
that is valid for a fully molecular gas, and 0.6 valid for 
a fully-ionized gas. For the opacities, we use the tables of 
Ferguson et al. (2005) for solar heavy-element abundances 
(set g7.02.tron, with the Grevesse & Noels abundances 
1993, X=0.7 and Z— 0.02). These opacities are displayed 
in Fig. IA.1I and compared to those of Bell & Lin (1994) 
generally used in these problems, which should be extrap- 
olated for low densities. 



Appendix B: Angular momentum provided by 
supernovae 

Rozyzka et al. assume that the disc receives all the linear 
momentum projected on the mid plane of the disc (i.e. 
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Fig. A.l. Ferguson et al. opacities used in this paper as a 
function of the temperature. The curves are labeled with 
the log of the density in g cm" 3 . For comparison the two 
thick lines show the Lin & Bell (dotted line) and the 
Ferguson et al. (solid line) opacities for the same density 
p=10~ 9 g cm" 3 . 



about Ptot/2). Thus the velocity reached by the remnant 
when it has expanded to a distance AR is 

Vc = ^ * (B.l) 

2 TrEAi? 2 V ; 

The remnant is stopped when this velocity equals 

the shear velocity at the edge relative to the explosion, 

flAR/2. Combining this value with Eq. IB. II gives Ai?fi na i 

and Wfinai, and finally the angular momentum given to the 

disc (see Rozyzka et al. for more details): 



A J = —PtotR. 
3n 



(B.2) 



From this value, Rozyzka et al. deduced that a rate of 
supernova explosions of 10" 4 to 10" 3 per year corresponds 
to the viscous transport in an a-disc having a=0.1 to 1. 

There is a big difference from CZ99. The Rozyczka et 
al. computation does not depend on the scale height, but 
only on the surface density which is given a priori. Taking 
the 3-dimensional aspect into account, CZ99 assume that 
the disk only receives the momentum corresponding to 
the velocity vectors inclined by less than H/R to the mid- 
plane. Calling 6* max the angle between the equatorial plane 
and the point where the remnant breaks through the disc 
surface, one gets the angular momentum supplied by a 
supernova, 



AJ = 



Ptot 
47T 



=i(0)2ircos(8)d0, 



(B.3) 
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For small angles, this expression can be approximated by 

(B.5) 



-^t-Sjinax 



where i? s ,max is the maximum extension of the supernova 
before it breaks vertically out of the disc: 



(3P tot \ 



1/4 



(B.6) 



The angular momentum given by Eq. lB.5l differs from that 
of Rozyczka et al. by a factor 9/4(77/i? s max ), which is gen- 
erally much smaller than unity; so our estimation is more 
conservative than that of Rozyczka et al. and it requires a 
larger number of supernovae to achieve the same momen- 
tum transport. 

The rate of supernovae per decade of radius is finally 
given by the relation: 



SN 



2M SN R 2 n 
AJlogioe 



(B.7) 



Appendix C: Equations of the clumpy disc 
structure 

We give here the complete set of equations used for the 
radial structure of the clumpy disc, though the majority 
of them have already been given in the text: 



^disp - (Gm d fi) 1/3 r- 

'cl 

m c i = lc\p 3 , 

where p is the density inside the clumps. 
M<jisc = 7T-R 2 E av . 
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where A is a factor of order unity, and Cav is the parameter 
( averaged on the disc: 



Cav — 

with 

Pav = 



47rG>a 

S av 
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dclV 

Ra J 



M = 37rS av d ' sp - v ';7 . (C.9) 

n (V) +1 



T = l c iK{p,T)pmax(l, l -^j. (CIO) 

rp4 

aT ° s = 3r/8 + 1/2 + l/4r (C ' U) 

where T e g — R~ 7 , with 7 imposed. These equations are 
completed by the two conditions: 

^ = given constantl (C12) 

^ = given constant2 (C.13) 
H 

or equivalently given coverage factor, and given l c \/d c \. We 
thus have a set of 13 equations with 13 unknowns, which 
can be solved by an iteration procedure. We also compute 
the Jeans temperature of the clumps, Tj oans , to check their 
gravitational stability: 

Tjcans= (10 2 V™ c1 Vp) 2/3 - ( C - 14 ) 
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Fig. 2. Dependence of several parameters on the radius, for the model corresponding to the quasar spectrum shown 
in Fig. [1] The parameters of the model are M = 10 9 M Q , Eddington ratio 0.5, a — 0.3, and a dependency of T e g 
versus R: T c g oc R~ 0AJ , with i?trans equal to 0.03 pc (300 i?schw, about i? mar g)- In the first figure (left and top panel), 
the optical thickness is displayed for all solutions, but the optically thin solutions have been suppressed in the other 
figures. The a-disc parameters are shown with big red dots, and the Q-disc with small blue dots. The junction between 
the a and the Q discs is shown by the big dark square. 
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Fig. 3. Same as Fig. [2j but for the model corresponding to the luminous Seyfert spectrum shown in Fig. Q] The 
parameters of the model are M = 10 8 M Q , Eddington ratio 0.1, a = 0.3, and a dependency of T e s versus R: 
T c ff oc R~ 0A7 , with i?trans equal to 0.003 pc (300 i?s c hw, about i? m ar g /3). Due to the additional heating, the a-disc is 
gravitationally stable up to 0.1 pc. 
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Fig. 4. Same as Fig. [2] but for the model corresponding to the weak Seyfert spectrum shown in Fig. [I] The parameters 
of the model are M = 10 7 M©, Eddington ratio 0.01, a = 0.3, and a dependency of T e g versus R: T e g oc R~ 0A7 , with 
-Rtrans equal to 0.003 pc (3000 i?schw, about i? marg /100). Due to the additional heating, the a-disc is gravitationally 
stable up to at least 0.006 pc. 
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Fig. 8. t V i SC and M as a function of the radius in the stationary (constant M) and in the non stationary cases for 
the 3 previous models. The black dots mark the position of the transition between the a-disc and the Q-disc in the 
stationary case. 



